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Abstract 

Field-tuned quantum tunneling in two single-molecule magnets coupled antiferromagnetically 
and formed a supramolecule dimer is studied. We obtain step- like magnetization curves by 
means of the numerically exact solution of the time-dependent Schrodinger equation. The steps 
in magnetization curves show the phenomenon of quantum resonant tunneling quantitatively. 
The effects of the sweeping rate of applied field is discussed. These results obtained from 
quantum dynamical evolution well agree with the recent experiment [W.Wernsdorfer et al. Nature 
416(2002)406]. 
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The macroscopic quantum phenomena in molecular magnets has become a very attractive 
researching field. Many properties of these nanometer-sized magnetic particles and clusters, 
such as Mn 12 (s = 10), Fe 8 (s = 10) and Mn 4 (s = |) systems, have been well studied M - ' 6 ' 
both experimentally and theoretically. Theoretically, studying the phenomenon of quan- 
tum resonant tunneling of these molecular magnets could be based on Landau-Zener (LZ) 
transitions' 7 ''^, or based on numerically the solution of the time-dependent Schrodinger 
equation Pl'I 10 !. In Landau-Zener model, the magnetization curves could be obtained in a 
static and approximate way. Recently a supermolecular dimer [Mn 4 ] 2 is reported to be 
synthesized successfully by Werndorfer et al.' 11 '. In this kind of supermolecular dimer, two 
single-molecule magnets Mn4 antiferromagnetic coupled each other, which results in its quan- 
tum behavior quite different from two individual Mn4 molecule without coupling. In this 
paper, we calculate magnetization curves of a supermolecular dimer [Mn 4 ] 2 ,a single- molecule 
magnets, by numerically exact solution of the time-dependent Schrodinger equation. 

Following Wernsdorfer et al., the model Hamiltonian of the supermolecular dimer [Mn 4 ] 2 

is 

H = H 1 + H 2 + JS 1 - S 2 , (1) 

where J is the weak antiferromagnetic supercharging coupling constant. H x and H 2 are 
Hamiltonian for two individual Mn 4 molecules in the supermolecular dimer. It is known 
that the model Hamiltonian of an individual Mn 4 molecule is 

Hi = -DSl + E(S 2 xt - S 2 yi ) - gfi B S zl h z (t),i = 1,2, (2) 

where D and E are the axial anisotropic constants. h z (t) is the applied sweeping field along 
easy axis. We can easily obtain the energy eigenvalues of whole Hamiltonian H (Figure 1). In 
experiment the sweeping rate of h z (t) is very slow, so we can simulate it as a step- increased 
field, which means the h z increases a value Ah z every r time step and keeps constant during 
the r time intervals. Note that we can not use a sweeping rate as slow as experiment due to 
the limitation of our computing time. However, we can obtain the key macroscopic quantum 
phenomena in our calculation with relatively high sweeping rate. In this paper, we select 
D = 0.72k, J = 0.1£;, [11] and E = 0.0317A: [6] . Dynamic evolution follows time-dependent 
schrodinger equation and can be calculated by 

|tf (f)) = |tf (t + nr)) = exp[-iH(t + (n - l)r) • r] \V (t + (n - l)r)> , (3) 
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Meanwhile, |^ (t)) can be expanded as 

S S 9 

= Yl Yl (*) i m i> m 2),^ = -, (4) 

mi=-sra2=-s 

where |m 1 ,m 2 ) are the eigenstates of H that is 

#o = - DS 2 zl + J^! • S z2 - 9f i B (S zl + S z2 )h z (t), (5) 

We assume the initial states to be at | — |, — |), and the whole evolution process can be 
obtained by equation (4) step by step. In Ref.' 11 ', Wernsdorfer report five points (Figure 4 
of Ref.) of resonant tunneling that result in the steps in hysteresis loops. They considered 
the first point is caused by the resonant transition from | — |, — |^ to| — |, |), and the fourth 
point is caused by those from | — |, — §)to|— §, §)■ However, under the model Hamiltonian of 
Equation(l) and Equation(2), the transitions of these points are quenched for a half integer 
spin due to the parity symmetry I 12 ' 'I 13 ' . Therefore, there must be some kind of transverse field 
components resulted from the influence of the environmental degrees' 6 ', such as hyperfine 
and dipolar couplings, and it can be approximated a Gaussian distribution with a width 
a = 0.035T for such additional transverse environmental field. In this paper, we simply 
assume it to be a constant and h x = 0.01T along x axis, but do not lose the essential 
physics, the macroscopic quantum phenomena. 

The magnetization along the z axis can be simply defined by M =< S zl + S z2 >■ In 
Figure 2, we plot the magnetization curve responding to a time-dependent applied field with 
a constant transverse field h x = 0.01T. There are three steps in the magnetization curve. 
In order to know the details of state transitions, the states (t))) near to two sides of 
resonant points are recorded in our simulation and they are shown in Table 1, where the 
occupied probabilities (|a mijjn2 (t)| 2 ) are neglected to zeros if they less than 1CT 3 . Therefore, 
we can get clear information of the process of evolution and transition. Figure 1 and Table 
1 show that the first step occurs at h z = 0.198T from | — |, — |) to | — |, |) (and||, — §)), and 
the second step occurs at h z = 0.731T from | — |, — |) to | — |, |) (and||,— |)). These two 
resonant points fit well to experimental results (i.e. the point 2 and the point 4 of Fig.4 in 
Ref.M). The third step occurs at h z = 0.812T from |-|,|> (and| §,-§)) to||, §) (or||,|». 
There is no step at h z ~ — 0.33T in our magnetization curve, but the experiment reports 
a point of resonant tunneling (the point 1 in Fig.4 of Ref.' 11 !). The reason is that we have 
used a too fast sweep rate in our simulation. We will interpret it more detail late. In our 
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magnetization curve, since the step at h z ss — 0.33T from | — |, — |) to| — §, §) does not occur, 
therefore the step at h z = 0.87T from |-|, §) to||, §) (point 4 in the Fig.4 of Ref. [nl ) can 
not occur naturally. 

In order to interpret why there is no step at h z = — 0.33T in our magnetization curve, 
we firstly consider a utmost-slow process. At any t value, the eigenstates and eigenvalues 
(Figure 1) of whole Hamiltonian H(t) can be calculated by 

tf(t)|$(£)> = £|$(£)>, (6) 

Note that no matter how weak it is, the system always interact with environment which 
cause dissipation. Therefore, in a very very slow process, we can assume that the state 
\^(t)) of system evaluating from a initial state |^(0)) will always relax to the ground state 
\®(E m m)) of H(t). Figure 3 shows the magnetization curve of this utmost-slow process. 
There are two steps at h z — 0.3363T and h z 0.3363T in the curve. It means that the 
point of resonant tunneling at h z xs — 0.3363T (point 1 of Fig.4 of Ref.' 11 ') occurs when 
the sweeping rate of applied field is very slow. In figure 2, the sweeping rate of applied 
field in our calculation is c = = jjps- = 1000 Tesla/s, which is much more larger than 
the ones in experiment (0.140 Tesla/s, 0.035 Tesla/s and 0.00ATesla/s) [11] . Due to the 
limitation of our computer time, we can not do the calculation for a sweeping rate of applied 
field as slow as the one in experiment. We now try to simulate the magnetization process 
(Figure 4) only in a very sharp range of h z with sweeping rates as slow as the ones used 
in experiment. In our figure 4, (a) is calculated over a range of h z from —0.3364 Tesla to 
—0.3362 Tesla with parameters r = 10~ 8 s and Ah z = 10~ 9 Tesla (the sweeping rate ~ 0.10 
Tesla/s); (b) is calculated over a range of h z from —0.336295 Tesla to —0.336275 Tesla 
with parameters r = 10 _8 s and Ah z = 10~ 10 Tesla ( the sweeping rate ~ O.OlTesla/ s); (c) 
is the combination of (a) and (b). A very clear step occurs at h z xs — 0.336283T point, and 
it shows that the slower sweeping rate induces the higher step occurred. The recorded states 
(Table 2) at transition point h z —0.3363 Tesla show that the resonant tunneling is from 
I - 1' — I) — §, |) (and||, — !))• All these results fit well with the results of experiment' 11 ' . 
Therefore, it clear show that the reason for no step at point about h z ~ — 0.3363Tes/a is 
from too fast sweeping rate of the applied field in theoretical simulation. There are some 
small oscillations in the magnetization curve. It is caused from quantum fluctuations. 

In conclusion, We have studied the phenomenon of quantum resonant tunneling in a 
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supermolecular dimer [Mn 4 ] 2 of single- molecule magnets by numerically exact solution of 
the time-dependent Schrodinger equation. We obtain step-like magnetization curves which 
demonstrate quantum tunneling quantitatively We have calculated and discussed the affect 
to steps caused by different sweeping rate of applied field. It shows that some steps can not 
occur at some resonant points when the sweeping rate is too fast, but they could appear when 
the sweeping rate becomes enough slow. At a very narrow region near resonant point, we 
slow down the sweeping rate of applied field, some quantum resonant tunnelling appeared 
in experiment can appear. Meanwhile, theoretical calculation show that more slow rate 
induces more higher transition step. The results of our calculation fit very well with the 
experiment' 11 '. Note that since we do not take into account the effects of dissipation caused 
by environment, the magnetization curves we obtain can not reach a reversal saturation 
value even if the applied field increase to infinitive value. Therefore, if we want to calculate 
a whole hysteresis loop, a proper mechanism of dissipation should be taken into account. 
This work is supported by National Natural Science Foundation of China. 
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caption 

Tablel: Occupied probabilities |a mijjrV2 (t)| 2 at spin states |mi,m 2 ) (Equation (5)) at 
some points of the evolution process (Figure 2). The values are neglected to zero if they are 
smaller than 0.001. 
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Table2: Occupied probabilities |a mijm2 (t)| 2 at spin states |mi,m 2 ) (Equation (5)) at 
some points of the evolution process (Figure 4). The values are neglected to zero if they are 
smaller than 0.0001. 
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FIG. 1: The 100 spin state energies of the model Hamiltonian (Equation (1)) as a function of 
longitude applied field. A weak transverse field h x = O.OlTesla is take into account. 
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FIG. 2: Magnetization relaxation of ground state (Magnetization curve reponse to a sweeping field. 
The sweeping rate is c = = = WOOTesla/s. 
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FIG. 3: Magnetization curve based on a utmost-slow process. It suppose that the state (^(i)) of 
system always relax to the ground state \^(E m i n )) of H(t). 
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FIG. 4: Magnetization curves response to slow sweeping fields over very sharp ranges of h z . The 
sweeping rates are the same order with that used in experiment (RefJ 11 ]). 
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